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Abstract. We present high-contrast imaging data on the disk of the classical T Tauri star TW Hya. The images were obtained 
through the polarimetric dilferential imaging technique with the adaptive optics system NACO. Our commissioning data show 
the presence of polarized disk emission between 0.1" and 1.4" from the star We derive the first Ks-band radial polarized 
intensity distribution. We show that the polarized intensity compares well to shorter wavelengths surface brightness observations 
and confirm the previously reported gradual slope change around 0.8". These results show the potential of the new polarimetric 
differential imaging technique at 8m-class telescopes to map the inner regions of protoplanetary disks. 
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1. Introduction 

The simulations of viscous disks with an embedded Jupiter- 
mass planet (see, e.g. D'Angelo et al. I2003I I predict dramatic 
changes in the disk structure: gaps are being opened and shock 
waves follow the path of the orbiting body. Detecting these sig- 
natures is an important step in confirming that our models for 
planet formation and disk evolution are right and would pro- 
vide a unique possibility to study a newly formed exoplanet 
and its interaction with the circumstellar material. 

However, the investigations of protoplanetary disks require 
high-contrast, high-spatial-resolution imaging very close to the 
bright central star. Although in recent years numerous groups 
studied disks around nearby stars at their planet-forming age 
with the HST, instrumental limitations do not allow to probe 
the disk structure at scales similar to our inner Solar System. 
The technique of polarimetric differential imaging (PDI) as de- 
scribed by Kuhn et al. J2001> is capable of tracing scattered 
(i.e. polarized) light from a circumstellar disk very close to the 
central star. This promising method enhances the contrast be- 
tween the disk and the star and, thereby, might be important to 
detect the detailed inner structure of protoplanetary disks. 

The basic idea of the PDI is to take the difference of two or- 
thogonally polarized, simultaeneously acquired images of the 



same object in order to remove all non-polarized light compo- 
nents. The non-polarized light consists mainly of the central 
star's light and the speckle noise, which is - close to the star - 
the dominant noise source of ground-based observations at op- 
tical and near-infrared wavelengths (Racine et al. ll999l l. After 
this subtraction only the polarized light, such as scattered light 
from the disk, is left over. 



To explore the potential of this technique for planet forma- 
tion studies as well as to demonstrate the NACO/VLT platform 
capabilities we targeted the 8 Myr old (Webb et al. I1999> and 
very close (56 + 7 pc, Wichmann et al. il998j classical T Tauri 
star TW Hya (RA: 11*'01'"51.9^ Dec:-34°42'17" J2000). At 
the visual wavelength regime TW Hya shows a polarization 
variability between 0-3% (Mekkaden il998» . Being an excel- 
lent target for planet formation studies, during the past few 
years several attempts have been made to image its disk as 
close as possible to the star and to identify possible disk struc- 
tures (Krist et al. 120001 Weinberger et al. I2002> . In addition, 
the presence of a giant planet was predicted (Calvet et al. 
1^021 ) from modeling the spectral energy distribution (SED) 
(see also Steinacker & Henning I2003I I. The first application 
of the PDI to TW Hya by Potter (200^ 200:^B made use of 
the 36 element curvature sensing Adaptive Optics (AO) system 
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* The data presented here were obtained during the commissioning 
of the NACO instrument at the Paranal Observatory of ESO. 



Hokupa'a (Graves et al. l2000l l. of the Gemini North telescope. 
This H-band data set probed the circumstellar material closer 
than 0.5"to the star. 



2 



D. Apai et al.: NACO Polarimetric Differential Imaging of TW Hya: 



In the following we describe the first results obtained with 
the PDI technique and the 188 element Shack-Hartmann AO 
system NACO, which is attached to the Nasmyth B port of the 
VLT UT4. 



2. Observations and Data Reduction 

The observations were carried out in April 2002 during the 
commissioning of the NACO/VLT system (Lenzen et al. . 19981 
Rousset et al. l2003l Hartung et al. l2000t in the Ks band. In 
order to simultaneously measure two orthogonal components 
of the polarized light, a Wollaston prism was introduced in the 
light path. The Wollaston prism splits the light into an ordinary 
and an extraordinary beam (o- and e-beam), separated by 3.5" 
in the Ks-band. To eliminate the instrumental polarization and 
to increase the signal-to-noise ratio (SNR) of the data, a redun- 
dant data set with Wollaston angles of 0°, 45°, 90° and 135° 
was acquired. Each of these images has a field of view (FOV) 
of 3" X 29.4" sampled with a pixel size of 0.027"/pixel. At 
every Wollaston position a 3-point dithering was applied with 
steps of about 9" along the field of view to allow sky subtrac- 
tion and to reduce the influence of bad pixels. No additional 
off-source sky frames were taken. 

To ensure the high dynamic range needed for disk analy- 
sis, we repeated the complete set of these observations with 
defined integration times of 0.4 s and 30 s. The total time spent 
on the source was 24 s and 1800 s in the short and long ex- 
posure series, respectively. As the TW Hya saturated the de- 
tector, the AO-performance can not be accurately measured. 
However, based on simulations with the NAOS Prepartion 
Software vl.62 we estimate a typical Strehl-ratio of 50% for 
these observations. 

The data reduction was carried out using self-developed 
IDL routines. First, every frame was manually inspected and 
those showing reflections or electronic ghosts were excluded 
from further reduction. Hot pixels were efficiently removed by 
a 3.5-sigma filtering process. The raw images were sky sub- 
tracted and flat field corrected. The sky frame was calculated 
individually for each group of images taken with the same po- 
larization angle. Following these basic corrections, the ordinary 
and extraordinary peaks of the individual exposures were cen- 
tered and extracted. To remove all of the non-polarized inten- 
sity (mainly from the central star) we subtracted images of or- 
dinary polarization from those of extraordinary polarization. 

The fine alignment of the individual images before subtrac- 
tion and co-addition was carried out by a two-level gaussian 
fitting procedure (see also Sectionl^. 

The subtraction process gives the orthogonal g, and U, 
Stokes components, ; standing for the number of dithering 
positions. For the 4 angles of the Wollaston prism we thus 
derived the values for Q'^, U'^^, Q^"' and C/.'^^ (see, e.g. 
Huard [T997b . Using the redundancy of our data set we derived 
the Q and U mean polarization vectors: 

1=1 1=1 



and the polarized intensity: 

This reduction procedure was repeated independently on 
both the long and the short exposure series, resulting in the 
^^shon and ^'^Long polarized intensity maps. 

Due to the brightness of TW Hya, the inner regions of the 
raw frames exceeded the regime of linear detector response. 
These data points (r < 0.06" in the short exposure series and 
r < 0.45" in the long exposure series) were excluded from the 
data reduction and further analysis. The field of view (FOV) 
was limited to ~ 3" by the field mask. During the dithering 
process slight movements perpendicular to the dithering direc- 
tions occured. Since the resulting FOV is the intersection of the 
FOVs of the individual frames, the final field of view was re- 
duced to ~ 0.54" and ~ 1.43" for the short and long exposure 
series, respectively. 

To characterize the surface brightness distribution of the 
disk we measured the radial profile of the polarized intensity, 
derived from 3-pixel-wide apertures, roughly equivalent to the 
full width at half maximum (FWHM) of the point spread func- 
tion (PSF). 

3. Signal Analysis 

For detecting the polarization signal as close as possible to the 
star a very careful data inspection is essential - residuals from 
the subtractions as well as sub-pixel misalignments between 
the individual frames are the dominant noise sources. While 
the scalar PI image contains no information on the direction 
of the polarization, the Q- and U-components are projections 
of the polarization vector and thus show a sinusoidal modula- 
tion with the angle in respect to the polarization direction of 
the Wollaston-prism (i.e. a maximum at the angles parallel to 
the polarization direction of the Wollaston prism and a mini- 
mum perpendicular to it). In the case of axisymmetrically dis- 
tributed dust this modulation leads to the butterfly pattern as 
shown in the right panel of Fig. ^ and described in Kuhn et al. 
J2001I I. By using the information about the polarization orien- 
tation and by comparing this pattern to what is expected from 
an axi-symmetric light scatterer, we can better discriminate be- 
tween noise and signal. 

The right panel of Fig. ^ shows the Q-component of the 
polarization signal radially integrated (over 5-degree bins) in 
an annulus between 0.75" and 1.0" from the long exposure se- 
ries. A very similar, but 7r/2-shifted modulation is evident in 
the U-images. The clearly sinusoidal modulation arises from 
the strong polarization signal. When probing the extent of the 
polarized circumstellar material, we looked for the presence of 
a similar pattern (phase- and frequency-correct sinusoidal mod- 
ulation) as a function of radial distance from TW Hya. 

We found an extended polarization pattern with a high 
signal-to-noise ratio between 0.5 and 1.4" distance from the 
star in both the Q and U long exposure images. In the short ex- 
posure images we identified light-scattering material between 
0.2" and 0.4" from TW Hya in both the Q- and U-components. 
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Fig. 1. Left paneLThe characteristic butterfly pattern of an extended polarized source around TW Hya in the Qhong image. This 
pattern can be identified between 0.5" and 1 .4" from the star on both the Q- and U (long exposure) images. The image is centered 
on TW Hya. Right panel: Counts in the Q-component of the polarization vector as a function of position angle in an annulus 
between 0.75" and 1" from the long exposure series. The asterisks mark the data points averaged over 5 degree bins, while 
the curve is the best fitting frequency-fixed cosine. The strong sinusoidal modulation indicates that the light-scattering dust is 
distributed nearly axisymmetrically around the light source. 




Fig. 2. Left panel: Same as left panel of Fig. 1 but for the the Qshon image. The polarization pattern can be identified between 0. 1 " 
and 0.4" on the short exposure image. Right panel: Same as Fig.^for an annulus between 0. 1" and 0.2" from the short exposure 
series. The sinusoidal modulation (with correct phase and frequency) shows the presence of light-scattering dust as close as 0. 1" 
to the star. For a comparison of the residuals from the PSF-comparison star plotted in the same scale, see right panel of Fig.|3] 



Inspection of the higher signal-to-noise Q-component (short 
exposure) showed this pattern also between 0.1"and 0.2" (see 
the left panel of Fig.|2j- 

In order to estimate the contrast ratio between the polarized 
and the unpolarized fluxes we compared the fluxes of the di- 
rect images (after sky subtraction and bad pixel removal) to the 



amplitude of the fitted cosine. The typical contrast ratios are 
Fpiu ^1% and Fpiu ^3.5-4% for radii of r = 0.1" - 0.3" and 
r ^ I", respectively. 
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Fig. 3. Left panel: TW Hya Short U image -pol. signal between 0.1" and 0.4". Right panel: Similar plot for the PSF comparison 
star GSC 07208-00319. No polarization signal visible (the plot range and analysis parameters are set to be the same as for the 
right panel in Fig.O. 



4. Reliability Tests and Confirmation 

To confirm our results, we analyzed the Q- and U-components 
of the PDl data set of the PSF-comparison star GSC 07208- 
00319, which was observed and reduced in an identical way to 
the TW Hya data set. In this case neither polarized emission in 
the PI image nor any sign of sinusoidal modulation in the Q- 
and U-band images could be found (see right panel of Fig.|3|l. 

We also conducted extensive tests to check the reliablity of 
our data reduction and analysis toolkit. During these test runs 
the effects of the following errors were investigated: sub-pixel 
shifts between the individual frames, different centering, inter- 
polation, and flat field normalization methods. For each test 
run the complete data reduction and analysis has been executed 
both for TW Hya and the PSF comparison star 

As the possible sub-pixel misalignments might produce a 
bipolar, positive-negative pattern mimicking the butterfly pat- 
tern of the extended polarization signal, special care was taken 
to investigate their effects on our results. In order to simulate 
centering problems, we have introduced normally distributed 
random shifts with FWHMs of 0.2, 0.5, 1.0 and 1.5 pixels after 
the centering step. The sinusoidal modulation of the TW Hya 
short exposure frames gets strongly distorted and asymmetric 
even for the sub-pixel misalignments (0.2 and 0.5) and the pat- 
tern is completely disrupted for misalignments as large as 1 or 
1.5 pixels. The short exposure images of the PSF comparison 
star start to show a sinusoidal modulation, but its frequency is 
one half of that of the polarization signal. These images show, 
in fact, offset negative and positive peaks which can easily be 
discriminated from the butterfly pattern shown in Fig. ^ The 
long exposure images are less sensitive to sub-pixel misalign- 
ment but are noticeably distorted in the case of large (1 or 1.5 
pixel) centering errors. 

In order to explore the influence of systematic errors in the 
centering procedure we have introduced fixed shifts of 0. 1 and 
0.2 pixels after our centering procedure. We found that sys- 



tematic centering errors as small as 0. 1 pixels lead to a strong 
positive-negative pattern and the disruption of the polarization 
signal. These systematic errors produced similar pattern for the 
PSF comparison star as reported for the case of random shifts 
and can, therefore, not resemble the correct polarization pat- 
tern. 

Additional tests included the comparison of centering 
methods based on cross-correlation and on Gauss fitting, bi- 
linear and cubic interpolations at the rotation of the images, 
as well as the local and global normalization of the flat field 
images. We found that the Gauss-fitting centering, the bilinear 
interpolation and the individual normalization of each stripe in 
the Wollaston-masked flat field leads to the highest signal-to- 
noise ratios. 

These test runs proved that the strongest noise sources in 
the reduction procedure are the subpixel misalignments and the 
flat field errors. We found no combination of the artificially in- 
troduced errors that could create a false (phase- and frequency- 
correct) sinusoidal modulation and thus could lead to the mis- 
interpretation of our data set. 

Chromatic effects from the Wollaston prism could also lead 
to slight differences between the o- and e-beams and thus to an 
imperfect subtraction. However, the MgF2 prism used in the 
CONICA camera introduces a relative displacement of only 
86 mas between the blue and red edge of the broad K-band 
(CONICA Manual). In the narrower Ks-band this effect is even 
smaller and is therefore not expected to influence the PDl ob- 
servations. 

As TW Hya and the PSF star were observed in very similar 
airmass ranges (1.02 - 1.08 and 1.09 - 1.22) and seeing con- 
ditions (0.6"), we expect that any atmospheric effect leading 
to a false polarization signal would effect the PSF comparison 
images in the same way. Since there is no evidence for such 
an atmospheric influence in the comparison images (see right 
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Radial Polarized Intensity Distribution 




Distance from TW Hya [ " ] 



Fig. 4. Relative radial polarized intensity distribution of the TW Hya disk between 0.5" and 1.4". The asterisks mark our 3- 
pixel-binned Ks-band PDI data, while the solid and dashed curves show the J- and H-band surface brightness distributions from 
Weinberger et al. J2002t . Although the surface brightness is not equivalent to the PI, applying an arbitrary offset shows the similar 
behaviour of the three curves. The slope change around 0.8" has been found already by Krist et al. ( I2000> and Weinberger et 
al. J2002t . The error bars represent the combined errors from the PSF-comparison residuals and the statistical errors of each 
annulus. 



panel of Fig. O, we exclude atmospheric chromatism as a poss- 
bility to mimick PDI signals. 

The examination of the data reduction, possible instrumen- 
tal and atmospheric effects led us to conclude that the extended 
polarization signal in our images indeed originate from the 
light scattering on the circumstellar dust around TW Hya. 

5. Results 

Based on the NACO observations, we obtained Ks-band im- 
ages of the circumstellar disk around the classical T Tauri star 
TW Hya, which probe the disk structure closer to the star than 
any previous observation. Our diffraction limited PDI of TW 
Hya shows an extended butterfly pattern (see Figs.mjji char- 
acteristic for a spatially resolved axisymmetric light scattering 
source. This pattern is present between 0.5" and 1.4" on the 
long exposure Q- and U-images and between 0.1" and 0.4" 
in the short exposure images. The observations of the regions 
closer than 0.1" to the star and between 0.4" and 0.5" do not 
have sufficient SNR for a reliable analysis. The SNR of the 



data set is sufficiently high to plot the azimuthally integrated 
Ks-band surface brightness density distribution between 0.5" 
and 1.4" (see Fig. EJ. 

6. Discussion 

6.1. Scattering 

The presence of the extended polarization signal across our 
FOV confirms the existence of light scattering material. The 
sinusoidal modulation of the Q- and U-images confirms the 
centrosymmetric polarization pattern. The fact that we detected 
such emission even as close as ~0.1" (~6 AU) to TW Hya is 
the first direct evidence that its dust disk extends so far in. This 
detection does not support the model of Krist et al. ( I2000t who 
proposed the presence of an inner 'dark' zone (r ^ 18 AU - 
0.3") with suppressed flaring angle, necessary to explain the 
exceptional brightness of the TW Hya disk. 

Based on the modeling of the SED, Calvet et al. J2002> 
predicted the presence of a mostly dust-evacuated inner gap at 
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r < 3 - 4 AU (r < 0.07"). Such a gap could be a tell-tale 
signpost of a massive planet in close orbit around TW Hya 
and would be of high importance. However, our commission- 
ing data do not allow to draw any definite conclusions on the 
existence of this feature. 

6.2. Radial Profile and Slope Change 

In Fig. 0] we show the comparison of the radial polarized 
intensity derived from our PDI observation and (arbitrarily 
shifted) radial surface brightness profiles measured by the HST 
(Weinberger 120021 . Although these measurements probe dif- 
ferent quantities at different wavelengths (Ks and J, H, respec- 
tively) a remarkable similarity is evident between the radial 
slopes. Based on this similarity, we conclude that the polar- 
ization degree is nearly independent of the radius and thus the 
polarized intensity characterizes the surface brightness. Thus, 
the PDI method can lead to contrast enhancement without the 
loss of information. 

The radial surface brightness profile between 0.9" and 1 .4" 
can be well fitted by a power-law function /ks,i °^ ^-3.i±o.3 
This is slightly steeper than the /j_h «: r"^-^*"-^ behaviour found 
by Weinberger et al. (t2002i) at 1.1 and 1.6 jum, for the region 
between 0.8" and 2.7". 

To place our results in context, we compare them to sur- 
face brightness slopes predicted by simple analytical approxi- 
mations (e.g. Whitney & Hartmann 19921. These approxima- 
tions consider only isotropic scattering and simple geometric 
effects and do not include wavelength-dependent terms. The 
predicted surface brightness power-law exponent for an opti- 
cally thick, geometrically thin disk is -3, while for a flaring 
disk an exponent of -2 is expected. Our Ks-band results in the 
0.9" < r < 1.4" region support the flat disk geometry or a 
very small flaring. However, the surface brightness slope is de- 
termined by an intricate interplay between many disk and dust 
parameters and is not a conclusive measure of the disk geome- 
try. 

As to underline this statement, a strong, gradual change in 
the brightness profile around 0.8" from TW Hya is evident. 
The annulus between 0.5" < r < 0.7" is characterized by the 
very gentle surface brightness behaviour of /ks,2 ~ ^-0 9±o.2 
The slope of the curve is even less steep than that expected 
from a flared disk (Whitney & Hartmann ll992t . 

Although the radial brightness slope is strongly influenced 
by the disk geometry, slope changes can have different origins. 
Even though a detailed analysis of the reasons of the slope 
changes is beyond the scope of the current paper, we refer to 
some of them. The most obvious one is an actual change in 
the surface density slope of the disk, translating into a slope 
change in the surface brightness. Such clear-cut changes might 
be the results of external or internal perturbation of the disk 
by a companion or a massive planet as seen in the case of HR 
4796A (see, e.g. Boccaletti et al. l2003> . 

A second cause might be the change of the flaring angle and 
through this the change in the illumination of the disk's surface. 
Alternatively, a radial dependence of the scattering cross sec- 
tion of the dust grains might introduce modulations in the sur- 



face brightness profile. Disentangling the effects of these pro- 
cesses requires high-angular resolution observations at longer 
wavelengths. 

We note, that this change in the slope at 0.8" has been de- 
tected previously at shorter wavelengths by HST direct imag- 
ing and coronographic studies (Krist et al. 120001 Weinberger et 
al. I2002t . together with other slope changes in the outer parts 
(>1.4") of the disk. 

Our data were acquired during the commissioning of the 
NACO instrument and suffer from high detector noise and lim- 
ited AO performance. The current performance of NACO is 
expected to result in roughly five times better SNR. 

The first results of a PDI survey of the TW Hya Association 
are presented by Huelamo & Brandner J2003t . 

For future PDI observers it might be important to point out 
some of the limiting factors in our data set. First, the dynamic 
range of the data set has to be increased by applying series of 
different exposure times. Furthermore, the flat field calibration 
plays a critical role and should be carried out with great care. 
Finally, the neccessarily saturated star will pose limitations to 
the accurate alignment of the individual frames, which will in 
practice determine the inner radius until which the polarization 
signal can be detected. 

7. Summary 

We presented the first high-contrast observations carried out 
with NACO/VLT using the PDI technique. The extended scat- 
tering pattern of TW Hya provides the first direct proof that the 
disk extends closer than 0.5", up to ~ 0.1" (~ 6 AU) from the 
star. 

We derive the first Ks-band radial polarized intensity profile 
between 0.5" and 1 .4". This profile strongly resembles the sur- 
face brightness profile seen at shorter wavelengths and shows 
as well a strong, gradual slope change around 0.8". The polar- 
ized intensity profile between 0.9" and 1.4" stands close to that 
expected from a flat disk. 

These results demonstrate the potential of the PDI tech- 
nique at the NACO/VLT platform for future planet formation 
studies. This technique is capable of imaging a new, before- 
hand unaccessible regime of circumstellar disks at very high 
resolution. Such images can help understanding the evolution 
of protoplanetary disks at the scales comparable to our inner 
Solar System. 
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